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The latest ALEPH measurements on heavy flavours are presented. In particular the measurements 
of \V„b\, the B(b — > S7) and a study on the width difference between mass eigenstates in the B s 
system are presented. 



I. INTRODUCTION 

The large sample of Z collected by ALEPH during the LEP1 consists of about 4.5 million hadronic decays. Out of 
them, about 1 million of bb events allow to investigate rare processes which could help constraining the CKM mixing 
matrix, before the start of the new experiments at the B factories and hadronic colliders. 

In this report I will present the measurements of the B(b — ■+ ulvX) jjj which allows a measurement of the \V u b\ 
CKM matrix element, the B(b — > sj) Q and a study on the width difference between mass eigenstates in the B s 
system. 

II. MEASUREMENT OF \V UB \ 

The measurement of \V u b\ is performed by measuring the semileptonic branching ratio of B hadrons into charmless 
final states (b — > ulvX ). In fact this rate is proportional to the V u f,\ 2 and can be computed in the framework of the 
Heavy Quark Expansion theory [gj . 

Previous measurements of \V u b\ were performed in both the exclusive and inclusive channels in the T(4S) [Q. 
The exclusive measurements show large theoretical uncertainties in the transition amplitude computation. Inclusive 
measurements overcame this difficulty looking for an excess of events at the endpoint of the lepton momentum 
distribution where the b — > c contribution vanishes. At LEP both methods are not viable, either due to the small 
branching ratio and large background in one case or to the small accuracy in reconstructing the B hadron rest frame 
in the other case. 

This measurement relies upon the fact that nearly 90% of the b — > ulvX decays are expected to have an invariant 
mass of the hadronic system Mx below the charm threshold of 1.87 GeV/c 2 . By contrast, only in 10% of the cases 
the lepton energy in the b rest frame, E* , is above the kinematic boundary for b — ► clvX transitions ||. 
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A. Event Selection 



Events are selected by requiring the presence of an identified lepton with momentum larger than 3 GeV/c. In the 
opposite hemisphere a b lifetime tag is applied to reduce non-b contamination to less than 2%. The neutrino energy 
and direction is estimated from the missing momentum in the lepton hemisphere with a typical accuracy of 280 mrad 
on its direction and 2 GeV on its energy. 

Particles coming from the charmless hadronic system X are reconstructed by means of two Neural Networks, one 
to select photons and the other charged particles. The B hadron rest frame is then reconstructed by adding the 
momenta of the lepton, the neutrino and the selected particles. The total energy is determined by assigning a mass 
of 5.38 GeV/c 2 to the total system. The momentum and angular resolution, obtained from the simulation, are 4.5 
GeV/c and 60 mrad respectively. 

The discrimination of the b — > ulvX signal events from the background b — > c transitions is made on a statistical 
basis and is based on the fact that the c quark is heavier with respect to the u quark, leading to different kinematical 
properties for the two final states. Because of resolution effects the separation that uses a single kinematical variable 
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such as Mx can be considerably improved by combining more information characterizing the leptonic and the hadronic 
parts into a multivariate Neural Network. 

The choice of the variables is based not only in such a way to have a good discrimination between signal and 
background, but also on a reduced sensitivity to the specific composition of the X u system. 

Among the Neural Network inputs, the ones with larger separating power are the invariant mass, the sphericity 
and the lepton momentum. The Neural Network output is close to 1 for signal events and close to for background 
events. 

The signal Monte Carlo use a Hybrid model following Ramirez et al. 0. At low hadronic energy (below 1.6 GeV), 
only resonant states are produced, while for large energy, non-resonant states are expected to contribute, which occurs 
in 75% of the cases. 

Figure |l| shows the Neural Network output for the Monte Carlo and data after all the selection cuts. The signal 
is extracted from a binned likelihood fit to the Neural Network output between 0.6 and 1.0, which gives the smallest 
total relative error on the measurement. 
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data and Monte Carlo with no b — > u transitions (points) compared to the b 
between the data and Monte Carlo with the fitted value of b — > u. 



u contribution (histogram), and d) difference 



The Monte Carlo is normalized to the same number of events in the data to reduce the sensitivity to the assumed 
efficiencies of the analysis cuts. The first bin is excluded in this normalisation procedure, to minimise the effects of 
uncertainties of the background events in the fit. In this way, the fit is effectively a fit to the shape of the signal and 
background distributions. 



B. Results and systematic checks 



An excess of 303 ± 88 events is clearly seen in the region where the signal is expected, and is in good agreement 
with the predicted Neural Network shape. 

From the fitted number of events one can extract the value of the branching ratio: 

B(b -> ulv) = (1.73 ± 0.55 stat ± 0.51 syst6 ^ c ± 0.21 syst6 ^ u ) x lO" 3 . (1) 
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The systematics come from the uncertainties in modelling the b — > c and the b — > ulvX transitions. For the b — ► c 
systematic source, the main uncertainties are the ones which come from the knowledge of the charm topological 
branching ratios and the statistical uncertainties on BR(b — > Z), BR(b — > c — > Z) and < Xb > (the average energy 
fraction of the B hadron) which have been measured by ALEPH |6) . The systematic coming from the modelling of 
the inclusive decays of the b — > u transitions dominates the second systematic error. 

Extensive checks have been made to support the measurement. The fit is stable against variations on the fit interval 
of the Neural Network output, as well as changing some variables in input of the Neural Network. 

Since neutral hadrons have not been considered when reconstructing the b hadron, a bad simulation of the b — ► c 
states involving energetic neutral hadrons would alter the background in the region of high values of the Neural 
Network output. The neutral hadronic energy distribution in a 30° cone around the lepton is different for final states 
with and without allowing a measurement of the inclusive production rate of K® in D meson decays. Good 
agreement is observed between data and simulation. The rate found is in very good agreement with the one measured 
by MARKIII. 

Finally, since no vertexing information has been used in the selection, a common vertex between the lepton and 
the charged hadronic system is searched for, with a cut on the \ 2 probability larger than 0.2. The efficiency for this 
cut is larger for b — > u transitions with respect to the b — > c, due to the presence of the charm vertex in the latter. 
This difference is further enhanced in the Neural Network region close to 1, due to the smaller charged multiplicity 
and poor vertex of the b — > c transitions. The good agreement of vertexing efficiencies in data and Monte Carlo is 
shown in Fig.||, where the effect of the signal inclusion is evident. 
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FIG. 2. Ratio of the vertexing efficiencies between data and Monte Carlo as a function of the Neural Net output cut. Black 
circles are with b — » u transitions, while open circles without. 



From this measurement and the average b lifetime one can extract the value of the \V u b\ matrix element, by using 
the relation obtained in the framework of the Heavy Quark Expansion Theory S: 



\V ub \ 2 = 20.98 g(6 7 n y' )L6pS (l ± 0.05 pert ± 0.06 m J x 10 
U.UUz tb 

where 77, = (1.554 ± 0.013) ps is the average b hadron lifetime: 

1Kb | 2 = (18.68 ± 5.94 stat ± 5.94 syst ± 1A5 HQE ) x lO" 6 . 
corresponding to \V ub \ = (4.16 ± 1.02) x 10~ 3 . 
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III. MEASUREMENT OF THE B(b -> sj) 



The b — > S7 decay is a flavour changing neutral current that in the Standard Model proceeds via an electromagnetic 
penguin diagram in which the photon is radiated from either the W or one of the quark lines and its branching ratio 
is predicted to be (3.76 ± 0.30) x 10~ 4 at the Z peak. 

Virtual particles in the loop may be replaced by non-Standard Model particles, such as charged Higgs bosons or 
supersymmetric particles which could either enhance or suppress the decay rate, making it sensitive to physics beyond 
the Standard Model 0. 

The b — y sj decay is expected to be dominated by two body decays, leading to the presence of a very energetic 
photon in the final state associated with a system of high momentum and rapidity hadrons, originating from a 
displaced secondary vertex. 

The composition of the signal b — > sj Monte Carlo is based primarily on predictions from the Heavy Quark Effective 
Theory as well as on the CLEO measurements of the inclusive b — > sj branching ratio and the exclusive B — ► K* (892)7 
§■ 

A. Event Selection 

The event reconstruction starts from the requirement of a photon with an energy larger than 10 GeV which does 
not gives any ir° when paired with other photons in the event. In that hemisphere, 7r° and K s mesons are searched for. 
The remaining tracks were assigned a probability to come from the b — > 57 according to their momentum, rapidity 
with respect to the B hadron direction and, in the case of charged tracks, their three dimensional impact parameter 
significance. The probability functions were derived from the simulation. 

The reconstruction of the b —y sj candidate is then made by the photon and the hadronic system in the same 
hemisphere resulting by adding the ir° and K s , charged tracks and neutral calorimetric objects in decreasing 
b — y S7 probability. The candidate is accepted if the jet mass lies within 700 MeV/c 2 of the mean B meson mass, the 
hadronic system has a mass smaller than 4 GeV/c 2 and a multiplicity smaller than eight objects. 

The hemisphere opposite to the candidate is required to be b-like using a lifetime b tag. 

After the preselection there remain 1560 hadronic events, which are then splitted into eight different categories 
depending on : 

• the value of the length of the major axis of the shower ellipse, 07, of the photon candidate in the electromagnetic 
calorimeter 

• the energy of the hard photon, E* 

• the b tag probability of the opposite hemisphere. 

The signal is extracted from a binned log-likelihood fit of the E* distributions of the eight subsamples, using the 
corresponding distributions for the signal and background simulations, taking into account the finite Monte Carlo 
statistics in each bin. The five parameters in the fit are iVj,-^, NpsRi Af(6— ►cW > -^(non-bW arm N ot her, which are, 
respectively, the total number of signal, Final State Radiation, (b — > c)ir°, (non — b)7T° and 'other' background events. 

Figure || shows the E* distribution for the purest signal sub-sample in data and Monte Carlo. The clear excess in 
data is nicely explained by the signal contribution. 

B. Results and systematic checks 

From the number of fitted events and correcting for the efficiency in the selection cuts the inclusive b — > S7 branching 
ratio is determined to be 
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B(b aj) = (3.11 ± 0.80 sta t ± 0.72 syst ) x 1(T 4 



(4) 



where the statistical error takes into account also the finite Monte Carlo statistics. The largest systematic error comes 
from the shape of the E* background distributions as well as on the relative proportion of each background source in 
each sub-sample. This is assessed by observing the change in the measured branching ratio as the boundaries between 
the eight sub-classes are varied. The uncertainty coming from the energy calibration has been estimated by varying 
both the ECAL and HCAL calibrations. Finally the uncertainty of the baryonic b — > 57 decays has been assessed by 
repeating the fit setting that fraction to zero. 
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FIG. 3. The energy of the photon in the rest frame of the reconstructed jet. The top figure shows the data and the Monte 
Carlo background events in the region of high purity signal. The bottom figure shows the excess in data after the subtraction 
of the Monte Carlo background and the signal distribution from the fit. 

Checks have been performed which are consistent with the expected b — > sj hypothesis: 

• there is evidence of lifetime in the same hemisphere as the photon 

• there is an excess of high momentum kaons 

• the shape of the excess 07 distribution is characteristic of single photons. 





IV. A STUDY ON THE WIDTH DIFFERENCE BETWEEN MASS EIGENSTATES IN THE B s SYSTEM 

One of the most challenging measurements in the B s — B s system is the measurement of the mass difference between 
the two mass eigenstates, AM S , which allows to constrain the CKM unitarity triangle. Since AM S is large || and 
therefore difficult to measure, a complementary insight might come from the width difference Ar s , which is related 
to AM S via the relation |[o| : 

AT S _ 3 m 2 b V Q cD , , 

s 2 m t Vqcd 

where the ratio of the QCD correction factors (rj) is expected to be of the order of unity, and does not depend on the 
CKM matrix elements |ll[]. The width difference is expected to be of the order of 10% of the total width. 

One of the simplest way to investigate the width difference is to measure directly one of the two components of the 
B s lifetime. The decay B s -> Di* )+ Di* } ~X is dominantly CP even [fL2| and, if one neglects CP violating effects, is 
the short lifetime eigenstate. 
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The investigation is performed in the case when each of the two D s decay into a <b. The <b are identified in the 
K + K~ decay mode. 



A. Event selection 

ALEPH has undergone a refinement of its tracking performances of the data already taken at LEP1. Firstly a 
better pattern recognition in the Silicon Vertex Detector allows to increase the vertexing efficiency of hadronic B 
decays of about 30%. In addition, thanks to the readout of the specific ionization induced by charged particles on the 
TPC pads, the dE/dx measurement efficiency has now increased from 85% to nearly 100%, without degrading very 
much the purity. This measurement is the first one which benefits from this upgrade. 

Each kaon candidate is required to have at least 1.5 GeV/c momentum and is identified by requiring the dE/dx to 
be consistent with the kaon hypothesis and vetoing pions. The angle between the two kaons in the <f> rest frame has 
to be larger than -0.8. The fab system is required to have at least 10 GeV/c in momentum and an invariant mass 
between 2.0 and 4.5 GeV/c 2 . Most of the combinatorial and fragmentation background is thus removed. 

Besides the signal events there remain other sources of double <b events: one or both <b can originate from fragmen- 
tation or combinatorial background. Three physics background are also present in Z — > bb and Z — > cc events: 

1. -B — > DD S (X), D s -» faX, D -f faX 

2. Br s ) — > D< S \X , £)( s ) — > faX and fa from fragmentation 

3. -D( s ) — > faX and fa from fragmentation 

The first reproduces exactly the signal signature and can not be removed in the selection. However the D decay is 
Cabibbo suppressed with respect to the signal process (D s -»• faX ~ 18%, D+ -> faX ~ 2.9%, D° ->• faX ~ 0.9%). 

To have a good tracking quality, each kaon must have at least one hit in the VDET. The two <b are then constrained 
to a common vertex with the \ 2 probability of the vertex greater than 1%. At this point it is possible to reconstruct 
the B s decay length as the distance of the fab vertex from the primary vertex projected along the momentum direction. 
To reject non b events a b-tag is also demanded. 

After this chain of cuts, the contribution due to <b from fragmentation is negligible. 

The global efficiency is about 8% with a b-purity of 84% in the final sample. 



B. Results 

The B s lifetime is determined from the proper decay time distribution of the fab events. For each B s candidate, the 
proper time is obtained from the decay length I of the fab system and the B s boost. The decay length is measured 
in three dimensions by projecting the vector joining the interaction point and the <p(b decay vertex onto the direction 
of flight of the fa 4 ) resultant momentum. The typical resolution of the fa 4 ) decay vertex along the direction of flight is 
200/im. The boost of the B s is computed from the nucleated jet method [[l3| starting from the two (b tracks. 

The B s is extracted from an unbinned likelihood fit to the proper time distribution of the B s candidates. The 
background events amount to 78% of the total, and their proper time distribution has been parametrised from the 
events in the two (b sidebands. 

Figure |^ shows the result of the fit: 

T short _ ( L42 ± o 23 ± 0.16)ps. (6) 
The main systematics comes from the combinatorial background shape parametrisation. 
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FIG. 4. Proper time distribution of the two <j> candidates for events selected in the data. 



From this measurement one could extract the difference in the width of the two B s mass eigenstates: 

Ar 1 - Tl hort 

= 2 -22 — = 0.24 ± 0.35 

r tb„ 



(7) 



where tb s is the average B s lifetime (1.61± 0.10 ps) and the statistical and systematic errors are combined. 
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